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Abstract. We interpret the three periods detected in OGLE LMC Cepheids SC3-360128 and SC5-338399 as 
corresponding to the first three overtones of radial pulsations. This interpretation imposes stringent constraints 
on parameters of the stars and on their evolutionary status, which could only be the first crossing of the instability 
strip. Evolutionary models reproducing measured periods exist only in a restricted range of metallicities {Z = 
0.004 — 0.007). The models impose an upper limit on the extent of overshooting from the convective core. Absolute 
magnitude of each star is confined to a narrow interval. This allows to derive a new estimate of the distance to the 
LMC. We obtain m — M ranging from 18™ 34 to 18™53, with a systematic difference between the two stars of about 
0^13. The rates of period change predicted by the models are formally in conflict with the derived observational 
limits, though the uncertainities of measured P may be underestimated. If the discrepancy is confirmed, it would 
constitute a significant challenge to the stellar evolution theory. 
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1. Introduction 

, Triple-mode radial pulsators are extremely rare. For many 
' years AC Andromedae, with the fundamental and the 
\ first two overtones simultaneously excited (Fitch & Szeidl 
' 1976), was the only known object of this type. Recently, 
two new triple-mode variables have been identified among 
short period Cepheids of the Large Magellanic Cloud 
(Moskahk, Kolaczkowski & Mizerski 2004, Moskalik & 
■ Kolaczkowski 2005). Both stars, LMC SC3-360128 (here- 
' after Star 1) and LMC SC5-338399 (hereafter Star 2), 
have previously been known to pulsate in two ra- 
dial modes, namely in the first and second overtones 
(Soszyiiski et al. 2000). The analysis of Moskalik & 
Kolaczkowski (2005) clearly shows presence of the third in- 
dependent periodicity. The period ratio of P3/P2 = 0.840, 
the same in both Cepheids, points towards identification 
of the newly detected mode as a third radial overtone. 
The two stars are the first known pulsators with three 
consecutive radial overtones simultaneously excited. 

Even if certain forms of pulsation are rare, it does not 
mean that they are unimportant. The best example is the 
role of double-mode Cepheid data in stimulating the revi- 
sion of stellar opacities. Each measured period of an iden- 
tified mode yields an accurate constraint on stellar pa- 
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rameters. For our objects we have three such constraints. 
Another advantage is the membership of both Cepheids in 
the LMC. The stars have been and will be systematically 
monitored. Moreover, we have reasonably good estimate 
of their absolute luminosities. This estimate is helpful but 
not precise, because the distance to the LMC is still de- 
bated. In fact, we will show that seismological analysis of 
the two triple-mode Cepheids yields an independent dis- 
tance determination. 



Our aim is to extract maximum information on the 
two stars from all the available data. After summarizing 
observational information for the two objects (Section 2), 
we focus on the inference from the values of measured 
periods. Our analysis in Section 3 is similar to that of 
Kovacs & Buchler (1994) for AC And. We conclude, that 
both triple-mode Cepheids must be in the evolutionary 
phase of the Hertzsprung gap crossing and we confront 
seismologically derived stellar parameters with predictions 
of the appropriate evolutionary models (Section 4). This 
further constrains the ranges of admissible stellar parame- 
ters. Section 5 is devoted to the LMC distance estimation 
implied by our models. In Section 6 we compare theo- 
retically predicted rates of period change with the limits 
inferred from observations. Finally, in Section 7 we sum- 
marize our conclusions and discuss future observations. 
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Table 1. Triple-mode Cepheids in the LMC 



OGLE No. 


Pi 


P2 


Pz 


V 


/ 


E{B-V) 


{V - /)o 


Wi 




[day] 


[day] 


[day] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


LMC SC3-360128 


0.541279 


0.436049 


0.366301 


17.53 


16.98 


0.134 


0.38 


16.13 


LMC SC5-338399 


0.579510 


0.466624 


0.392116 


17.37 


16.87 


0.133 


0.33 


16.09 



2. Data 

Main observational properties of both triple-mode 
Cepheids are summarized in Table 1. Pulsation periods 
are determined from OGLE J-band photometry reduced 
with the DIA algorithm (Zebrun ct al. 2001). All peri- 
ods are accurate to better than 7 x 10~^day. Their er- 
rors are not quoted here; we believe they are smaller than 
uncertainities of the model periods and, therefore, they 
are inconsequential for our analysis. The mean V^-band 
and 7-band magnitudes are taken from Soszynski et al. 
(2000). The average E{B - V) reddening for the field of 
each Cephcid (Udalski et al. 1999) is used to obtain dere- 
denned {V — I)o colours. Standard extinction coefficients 
are adopted: Ay = 3.24 E{B - V),Ai = 1.96 E{B - V). 
In the last column of the table we give Wessenheit index, 
W/ = 7—1.55 X (V — I), which is an extinction insensitive 
brightness indicator (Madore & Preedman 1991). 

3. Constraints on stellar parameters from 
pulsation periods 

The basic assumption behind our analysis is that the three 
periods detected in the triple-mode Cepheids correspond 
to the first three consecutive overtones of radial pulsa- 
tions. The constraints on the stellar parmeters are derived 
by equating the three measured periods Pn.o with periods 
calculated in plausible stellar models, Pn,c- In this section 
we use only envelope models, extending from the photo- 
sphere down to the temperature of logT = 7.2. This is 
the standard approach in studies of radial pulsations in 
giants. We checked that the calculated periods are suffi- 
cently insensitive to the location of the inner boundary of 
the envelope. The model periods depend on six parametrs. 
These are: the mass, M, luminosity, L, effective tempera- 
ture, Teff, the fractional abundances of hydrogen, X , and 
of elements heavier than helium, Z, and the mixing length 
parameter, aeon- The pulsation data yield for each star 
three equations of the form 

Pn,c{M, L, TeS, X, Z, Q!con) = Pn,o 

and this is not enough to determine the values of all the 
unknown parameters. In principle, we could use mean 
magnitudes and colours of the stars, but these are un- 
certain. Therefore, we decided to find solutions for M, L 
and Teff for selected X, Z, and aeon in admissible ranges, 
relying only on the measured periods and requiring that 
all throe modes are linearly unstable at specified param- 
eters. Comparison with the observed mean colours and 
magnitudes is done a posteriori. 



Our static stellar models, both envelope models used 
here and the evolutionary models used in Section 4.2, are 
calculated with a modernized version of the Paczyriski 
(1970) codes. The main improvements consist in the im- 
plementation of the OPAL cquation-of-statc and opacity 
data (Rogers, Swenson & Iglcsias 1996; Iglesias & Rogers 
1996). Pulsation periods arc calculated with a modified 
version of the Dziembowski (1977) code. 

There are uncertainties in the calculated periods. We 
rely on the linear nonadiabatic approximation. We believe 
this is the most severe simplification. Comparison of pe- 
riods derived from linear and from nonlinear codes has 
been published only for monoperiodic RR Lyr models (e.g. 
Kovacs & Buchler 1988) and long period Cephcid models 
(e.g. MoskaUk & Buchler 1991; Fokin 1994). It is not clear 
how to scale these results to higher overtones in multiperi- 
odic ^ 3A'/0 Cepheids. 

Errors may result from using the mixing length the- 
ory to calculate convective flux in our static models and 
ignoring its Lagrangian perturbation in pulsating models. 
We believe that this is less essential and that our calcu- 
lations made with aeon = and aeon = 1 reflect well the 
uncertainty. The use of the Edddington approximation in 
the optically shallow layers is also an inadequacy but of 
secondary importance for the present application. 

Simultaneous fit to three radial mode periods can 
be achieved for several different sets of parameters {e.g. 
AC And; Kovacs & Buchler 1994). In order to assure that 
we find all the solutions we apply the following procedure: 

1. first, we select the chemical composition {X, Z) and 
the mixing length aeon- Keeping these three parame- 
ters fixed, we calculate a grid of stellar models cover- 
ing the mass range of 2 < M/Mq < 5 and luminosity 
range of 2 < log(L/Lo) < 3. 

2. at every point of the (M, L) grid we find a model which 
fits the observed period of the first overtone, Pi^o- This 
can be easily achieved by adjusting the last remaining 
parameter of the model, namely its effective tempera- 
ture. The solution for Tes is unique. The model fitting 
Pi docs not necesarily reproduce the other two peri- 
ods. We calculate period mismatches AP2 = P2,c—P2,o 

and AP3 = P3,e - P3,o. 

3. we find on the (M, L) plane the lines of APj = 0. 
Models located on these lines reproduce both Pi and 
P2. Similarly, we find the lines of AP3 — 0. The values 
of (M, L) at which the two types of lines cross each 
other define the models which simultaneously repro- 
duce all three periods. 
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Fig. 1. Simultaneous fit to three pulsation periods observed in Star 2. Three panels display models calculated for Qcon ~ 0, 
X = 0.70 and metaUicities of Z = 0.009, Z = 0.007 and Z = 0.005. Filled and open circles denote Unes of AP2 = and 
AP3 = 0, respectively. Different triple-mode solutions are marked A, B and C. Dotted area is the domain where all three modes 
are linearly unstable. 

Table 2. Seismic envelope models for Star 1 



Q^con 


X 


Z 


M/Mq 


log(L/L0) 


logTeff 


iV-I) 


m- M 












[K] 


[mag] 


[mag] 


Q 


0.74 


0.006 


3.29 


2.311 


3.814 


0.50 


18.43 








3.63 


2.403 


3.829 


0.43 


18.52 





0.74 


0.005 


3.15 


2.358 


3.828 


0.44 


18.42 








(3.67 


2.487 


3.849) 






1 


0.74 


0.007 


3.48 


2.299 


3.807 


0.52 


18.47 








3.67 


2.359 


3.818 


0.48 


18.52 


1 


0.74 


0.006 


3.21 


2.316 


3.817 


0.48 


18.42 








3.68 


2.423 


3.834 


0.41 


18.53 


1 


0.74 


0.005 


3.14 


2.376 


3.833 


0.42 


18.41 








3.59 


2.461 


3.845 


0.37 


18.51 


1 


0.74 


0.004 


3.14 


2.433 


3.847 


0.36 


18.42 








(3.60 


2.547 


3.865) 









0.70 


0.008 


3.28 


2.195 


3.787 


0.60 


18.39 








3.48 


2.263 


3.800 


0.55 


18.45 





0.70 


0.007 


3.10 


2.215 


3.795 


0.57 


18.37 








3.63 


2.366 


3.822 


0.46 


18.50 





0.70 


0.006 


2.98 


2.279 


3.813 


0.50 


18.36 








3.66 


2.437 


3.839 


0.39 


18.52 





0.70 


0.005 


2.88 


2.346 


3.832 


0.42 


18.35 








(3.66 


2.506 


3.856) 






1 


0.70 


0.007 


3.21 


2.240 


3.799 


0.56 


18.40 








3.67 


2.388 


3.827 


0.44 


18.51 


1 


0.70 


0.006 


2.96 


2.300 


3.819 


0.47 


18.36 








3.65 


2.440 


3.840 


0.39 


18.51 


1 


0.70 


0.005 


2.90 


2.362 


3.836 


0.40 


18.35 








(3.61 


2.487 


3.852) 






1 


0.70 


0.004 


2.86 


2.410 


3.848 


0.35 


18.35 








(3.63 


2.575 


3.873) 
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Table 3. Seismic envelope models for Star 2 





X 


Z 


M/Mq 


log(L/L0) 


logTeft 


(V-/) 


rn - M 












|K] 


[mag] 


[magj 





0.70 


0.008 


3.45 


2.333 


3.807 


0.52 


18.52 








3.65 


2.382 


3.815 


0.49 


18.57 





0.70 


0.007 


3.29 


2.367 


3.818 


0.48 


18.50 








O '7/1 
6.(4: 


Z.4Di 


6.000 


U.4z 


1 o art 
18. bU 





0.70 


0.006 


3.21 


2.424 


3.834 


0.41 


18.49 










Z.OZO 


Q o /I o^ 
0.848] 






1 


0.70 


0.008 


3.41 


2.337 


3.809 


0.51 


18.51 








3.74 


2.416 


3.823 


0.46 


18.58 


1 


0.70 


0.007 


3.27 


2.385 


3.824 


0.45 


18.49 








O. (0 




O.OOD 


n Art 
U.4U 


18. 09 


1 


0.70 


0.006 


3.25 


2.438 


3.837 


0.40 


18.50 








3.67 


2.504 


3.845 


0.36 


18.59 


1 
1 


u. / u 


u.uuo 


O.ZD 


9 /ISO 


Q 8/10 


U.oO 


lo.Oi 








(3.64 


2.563 


3.860) 









0.72 


0.007 


3.52 


2.397 


3.822 


0.46 


18.54 








3.66 


2.430 


3.827 


n A A 
0.44 


18.58 





0.72 


0.006 


3.38 


2.438 


3.834 


0.41 


18.53 








(3.72 


2.508 


3.845) 






1 


0.72 


0.007 


3.42 


2.399 


3.824 


0.45 


18.53 








3.73 


2.456 


3.832 


0.42 


18.60 


1 


0.72 


0.006 


3.45 


2.458 


3.838 


0.39 


18.54 








3.59 


2.482 


3.841 


0.38 


18.57 



In Fig. 1 we display the results of our procedure applied 
to Star 2, for aeon = 0, X = 0.70 and three different 
metahicity values of Z = 0.009, 0.007 and 0.005. The filled 
and open circles define the lines of AP2 — and AP3 = 0, 
respectively. The dotted area in each plot corresponds to 
models with all three modes linearly unstable. We recall 
here, that in the absence of resonances (as in our case), 
the linear excitation of all modes is the necessary condition 
for the triple- mode pulsations to exist. Among all models 
matching the observed periods, only those located in the 
dotted area correspond to acceptable solutions. 

As we see from Fig. 1, the problem of matching periods 
of a triple-mode Cephcid can have up to three different so- 
lutions. We call them A, B and C. Solution A exists for 
all choices of X, Z and aeon, but all three modes arc al- 
ways linearly stable. We will not discuss this solution any 
further. Solutions B and C can exist only for metallicities 
below a certain critical value. This is clearly seen in Fig. 1: 
for Z = 0.009 the hues AP2 = and AP3 = are close 
to each other, but do not cross. As the value of Z is low- 
ered, the two lines come closer and finally cross each other, 
generating solutions B and C. For the case displayed in 
Fig. 1, this happens for Z 0.0082. When Z is lowered 
further, first solution C and then solution B migrates out 
of the doted area. This sets a lower limit for the metallic- 
ity. Below this limit solutions B and C do exist, but do 
not correspond to triple-mode pulsations, because at least 



one of the modes is linearly stable. For the case of Fig. 1 
this lower limit is at Z 0.0056. 

We have performed search for solutions for both triple- 
mode Cepheids, exploring the range of metallicities from 
Z = 0.004 to 0.010 (in some cases to 0.020) and of hydro- 
gen abundances from X — 0.70 to 0.74. For each chemical 
composition we have calculated models with Ofcon = and 
with aeon = 1- The results of this survey arc summarized 
in Table 2 for Star 1 and in Tabic 3 for Star 2. Wc list only 
those combinations of X, Z and aeon for which at least 
one acceptable sohition exists. For low metallicities, solu- 
tion C may already fall outside the "all- modes-unstable" 
domain, in such a case it is given in parenthesis. For all ac- 
ceptable solutions wc provide the values of (V — I) colour 
and the distance modulus to the star, m — M . The latter 
is obtained by comparing Wessenheit index, Wj, derived 
from the models with the mean observed values, as given 
in Table 1. The conversion from model parameters to V- 
band and /-band magnitudes is done with tabular data 
based on Kurucz's (1998) stellar atmosphere models. 

The triple-mode pulsations strongly constrain parame- 
ters of the studied Cepheids. We find, that metallicities of 
both stars have to be in the range of Z = 0.0035 — 0.008, 
with limiting values somewhat dependent on X and as- 
sumed aeon- For Star 1, fitting the observed periods is 
possible for all values of X, but for Star 2 hydrogen abun- 
dance of X < 0.72 is required. The inferred stellar masses 
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are limited to ranges of 2.8 < M/Mq < 3.7 for Star 1 
and 3.2 < M/Mq < 3.8 for Star 2. The respective dis- 
tance moduli are in the ranges of 18™35 — 18™53 and 
18™49 - 18^60. 

In the next section we will show that the constraints 
imposed by triple-mode pulsations can be considerably 
tightened by combining the seismic envelope models with 
the appropriate evolutionary tracks. 
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4. Evolutionary models 

4.1. Evolutionary status of the stars 

The very short periods and low luminosities of the triple- 
mode Cepheids cannot be explained with the current 
models of stars in the phase of central helium burning. 
Whether or not a star enters the instability strip during 
this evolutionary phase depends on the extent of the blue 
loop. The latter is mainly determined by the mass and the 
metallicity of the star. Alibert et al. (1999) find that at 
Z = 0.004 their models with masses of 3.0 — 3.5Mq do en- 
ter the instability strip. However, the minimum value of 
Pi, which occurs for M = 3Mq is about 0^9. At Z = 0.01 
they find models with blue loop extending to the instabil- 
ity strip only at M > AMq. Such models have, naturally, 
much higher luminosities and pulsation periods than our 
stars. Alibert et al. (1999) assumed negligible overshoot- 
ing in their calculations. 

The effect of finite overshooting was taken into account 
in an extensive set of models of Girardi et al. (2000). In 
this case, the parameters we have derived for the triple- 
mode Cepheids are even more difficult to reconcile with 
the models in the helium burning phase. At Z = 0.004 
and M = 3Mq the loop does not extend to the instabil- 
ity strip. At M = 3.5Af0 the loop ends within the strip, 
but then the minimum value of Pi is about 1^15. Again, 
increasing Z only worsens the discrepancy. 

Determination of the extent of the blue loop is a subtle 
problem. Here we find the largest discrepancies between 
various authors. A more robust result is the difference of 
A log L ~ 0.2 between the first and the second crossing of 
the instability strip by stars in the 3 — AMq mass range. 
The logL values we find for the triple-mode Cepheids in- 
dicate that these two objects are caught during the short 
phase of the first crossing of the strip. 

4.2. Models in the first crossing of the instability strip 

We have computed a number of evolutionary tracks cover- 
ing the Hertzsprung gap crossing phase. From each track, 
we may determine L in the moment when the first over- 
tone period is equal to Pi^o- The value of L is virtually 
independent of aeon, but depends strongly on ctov, which 
is the ratio of the overshooting distance from the convec- 
tive core to the pressure distance scale calculated at the 
edge of the core. From our model calculations we find 



2.6 7 




0.48 0.5 0.52 0.54 0.56 

log (M/M„) 

Fig. 2. Evolutionary M ~ L relation for models with Pi = 
0'?541. Lines are derived for stellar models undergoing first 
crossing of the instability strip. Models are calculated assum- 
ing no overshooting from the convective core (qov = 0) during 
the main sequence evolution. Symbols represent seismic so- 
lutions for Star 1, derived with (non-evolutionary) envelope 
models (see Table 2). Filled and open symbols correspond 
to Ocon ~ 1 and 0, respectively. Pentagons correspond to 
X = 0.70, Z = 0.006, circles to X = 0.74, Z = 0.006, triangles 
to X = 0.74, Z = 0.005, and square to X = 0.74, Z = 0.007 
(there is no solution for aeon = in the last case). 



At specified values of A/, X , Z and Pi, the minimum value 
of L corresponds to aov = 

In Fig. 2 we display the M — L relations for no over- 
shooting evolutionary models of Star 1 (Pi — 0'?541). 
Different lines correspond to various chemical composi- 
tions. In the same plot we also display parameters of sev- 
eral seismic models of Star 1, also for different values of X 
and Z (see Table 2). We plot only models corresponding 
to solution P, because only for this solution agreement 
between seismic and evolutionary models can be reached. 
The lower panel of the plot explores the effect of varying 
stellar metallicity. For Z = 0.005 the seismic models do 
not agree with aov — evolutionary tracks. The agree- 
ment may be easily reached if overshooting of Ofov '^-^ 0.1 is 
allowed, which shifts the evolutionary M — L relation up. 
For Z = 0.006 the seismic models and the evolutionary 
tracks are almost in agreement. For Z = 0.007 the seis- 
mic model is placed below the corresponding evolutionary 
M ^ L line and no agreement between the two is possible. 
Fig. 2 illustrates the general pattern of model behaviour: 
with X, aeon and aov fixed, the exact match between seis- 
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Table 4. Triple-mode Cepheid models consistent with no overshooting evolutionary tracks 



X 


Ocon 


z 


M/Mq 


log(L/L0) 


logTeff 
[K] 


iv-i) 

[mag] 


m - M 
[mag] 


P/P 
[Myr-i] 


Star 1 


0.70 



1 


0.0061 
0.0063 


2.99 
3.02 


2.273 
2.280 


3.811 
3.813 


0.51 
0.50 


18.36 
18.37 


1.51 
1.53 


0.74 




1 


0.0057 
0.0062 


3.25 
3.26 


2.323 
2.308 


3.818 
3.814 


0.48 
0.50 


18.43 
18.43 


1.65 
1.61 


Star 2 


0.70 




1 


0.0065 
0.0069 


3.24 
3.26 


2.394 
2.392 


3.826 
3.826 


0.44 
0.44 


18.49 
18.49 


2.15 
2.18 


0.72 





0.0061 


3.39 


2.432 


3.832 


0.42 


18.53 


2.11 




1 

















mic models and evolutionary tracks can be reach only at 
one specific metallicity value. The matching values of Z 
becomes lower, when the overshooting distance is in- 
creased. 

The upper panel of Fig. 2 shows that match between 
seismic and evolutionary models is possible for any rea- 
sonable value of X, but yielding different values of stellar 
mass: M - ?,.2Mq for X = 0.74 and M ~ S.OMq for 
X = 0.70. Interestingly, metallicities of matching models 
are almost the same independently of X. 

In Table 4 we list parameters of the two stars resulting 
from matching seismic models to the evolutionary tracks 
calculated for aov = 0. The values obtained with two 
choices of aeon give some assessment of uncertainity due to 
poor knowledge of convection. We regard the parameters 
derived at higher X to be more likely, because we do not 
expect any significant Ho enrichment for LMC stars with 
M ~ 2>Mq. In the last two columns we provide the values 
of {y — I) colour and of the distance modulus to the stars, 
obtained in same way as in Tables 2 and 3. Bearing in mind 
uncertainty in the mean deredenned colour, which accord- 
ing to Udalski (private communication) is about 0™1, we 
regard the agrecnicmt in {V — I) to be quite satisfactory. 
The agreement may be even better, if smaller LMC reden- 
nings of Subramaniam (2004) are used. We will discuss the 
inferred values of the distance modulus later in this paper. 
In the last column of Table 4 we list the calculated rates 
of period change. These numbers will be discussed later, 
too. 

We stress, that with exception of P/P, all the values 
listed in Table 4 are derived from the envelope models. The 
requirement of matching the seismic envelope models and 
the evolutionary tracks was used only to a) reject solution 
C and b) to find the value of Z at which the exact match 
is achieved. 

4.3. Constraints on the extent of overshooting 

At fixed X and aeon, there is a minimum metallicity 
value allowed by seismic models. This lower limit for Z 



corresponds to an upper limit for the overshooting dis- 
tance, aov In Table 5 we give the maximum overshoot- 
ing distances for both triple-mode Cepheids. We obtain 
a stronger contraint for Star 2, especially if the higher 
value of X is adopted. Despite very similar masses, we 
cannot argue that the two stars should have similar 
They may well differ in the rotation rate and rotation in- 
duced mixing can mimic the convective overshooting. In 
any case, the limit of aov < 0.12 obtained for Star 2 is 
very interesting. 

The triple-mode Cepheids can be quite important for 
constraining element mixing in moderate mass stars. To 
demonstrate this point, we list in Table 5 also the stel- 
lar parameters of the two Cepheids derievcd at the ex- 
treme values of aov. These can be compared with param- 
eters derived for aov = (Table 4). As overshooting is in- 
creased, the matching models become less metal abundant 
and bluer. Thus, the limits on the extent of overshooting 
can be further improved once we precisely measure either 
metallicities or the mean deredened colours of the two 
stars. We assume here, that the composition in the enve- 
lope and in the photosphere is the same, which is valid 
unless the star is a fast rotator. 

5. The distance to LMC 

The values of distance moduli listed in Tables 4 and 5 are 
all between 18™34 and 18™53, i.e. well within currently 
considered range for the LMC. They are indeed very close 
to estimates based on analysis of RRd stars, which yields 
18^53 according to Kovacs (2000b) and 18™42 according 
to Popielski et al. (2000). From the double- mode Cepheid 
data, Kovacs (2000a) obtained the value of 18^54. All 
these determinations rely on fitting the observed pulsation 
periods, but for different objects and in different ways. 

For the same model parameters {X, Ocon), there is 
a systematic difference of about 0^13 between distance 
moduli derived for Star 1 and Star 2. We believe, that 
large part of this difference can be attributed to Wj 
measurement errors. Indeed, judging from faint Cepheids 
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X 


aeon 


z 


OCov 


M/Mq 


log(L/L0) 


logTeff 


iV-I) 


m - M 


P/P 














[K] 


[mag] 


[mag] 


[Myr-i] 


Star 1 


0.70 





0.0043 


0.33 


2.81 


2.393 


3.845 


0.36 


18.34 


1.69 




1 


0.0036 


0.33 


2.84 


2.428 


3.853 


0.33 


18.34 


1.81 


0.74 





0.0045 


0.22 


3.07 


2.392 


3.838 


0.40 


18.40 


1.61 




1 


0.0038 


0.19 


3.14 


2.443 


3.849 


0.35 


18.42 


1.75 


Star 2 


0.70 





0.0056 


0.12 


3.18 


2.450 


3.841 


0.38 


18.49 


2.07 




1 


0.0049 


0.07 


3.26 


2.493 


3.850 


0.34 


18.51 


2.26 


0.72 





0.0058 


0.05 


3.35 


2.452 


3.838 


0.39 


18.52 


2.12 




1 



















{Wi > 15™0) with independent photometry from two 
OGLE fields, the error of the mean Wessenheit index 
is about CP06 on average, but sometimes can exceed 
0^09. We cannot exclude, however, that some part of the 
CK"13 difference might be real. As shown by Subramaniam 
(2003), the m — M displays considerable variation even 
within the LMC bar. 

We want to sticss, that distances of the triple- mode 
Cepheids can be estimated very accurately once photo- 
metric precision is improved. As seen from Tables 4 and 

5, the total range of m — M allowed by seismic models is 
only 0^09 for Star 1 and 0^"04 for Star 2. In principle, the 
individual distances of the two stars can be determined 
with, respectively, ±2% and ±1% accuracy. 

6. Period changes 

The first crossing of the instability strip is by two orders of 
magnitude faster than either the second or the third cross- 
ing. Measuring the rates of period change should, there- 
fore, provide a crucial test for the triple-mode Cepheid 
models. 

The rates of period change, P/P, inferred from our 
models arc given in the last column of Tables 4 and 5. They 
are about 2Myr~^, somewhat below this value for Star 1 
and somewhat above it for Star 2. The rates are calculated 
for the first overtone, but they arc approximately the same 
for all modes. With such values of P/P, the (O — C) should 
reach 6 — 10% of the pulsation period in ten years. 

We do not have yet ten years worth of data at our 
disposal, but we will try to measure P with photometry 
already at hand. Our method consists in the least square 
fitting of the photometric measurements with the Fourier 
sum of the form 

I{t) = (J) -h V Ai cos[(u;i + -LOit) t + (l)i] + 

+ harmonics + combination terms 



where Wj = l-n/Pi are pulsation frequencies and uji = 
—2-KPi/Pf are the rates of frequency change (see e.g. 
Kepler 1993 for derivation). The nonlinear least square fit 
is performed with the Marquardt algorithm (Press et al. 
1986). In the following, we will discus only the rates of 
period change for the first overtone. The other two modes 
have much lower amplitudes and consequently their P is 
determined with 3-10 times larger errors. 

The method, when applied to DIA-reduced OGLE I- 
band photometry yields Pi/Pi = -1-1.7 ± 3.7Myr-i for 
Star 1 and A/A = +5.3 ± 6.4 Myr^^ for Star 2. Clearly, 
the OGLE timebase of ~ 1200 day is not sufficient to ob- 
tain significant P measurement. 

A considerably longer timebase (~ 2700 day) is ofered 
by MACHO photometry (Allsman & Axelrod 2001). Both 
triple-modes Cepheids have been observed by MACHO 
Project (LMC SC3-360128 = MACHO 77.8029.310; LMC 
SC5-338399 = MACHO 80.7074.120) and over 1500 mea- 
surements have been accumulated for each of them. 
Unfortunately, MACHO data arc of much lower quality. 
We have decided to perform initial filtering of MACHO 
photometry, by rejecting a) all outlying measurements 
deviating by more than 5cr from average brightness of 
the star and b) all measurements with photometric er- 
rors (given by MACHO) larger that 5 times the averge 
for the dataset. Further data clipping is done during the 
fitting procedure, when all points deviating by more than 
5cr from the fit are also rejected. For the final determina- 
tion of P we use the MACHO B-h&nd photometry. The 
analysis yields 

Pi/ Pi = -1.66 ± 0.83 Myr"^ for Star 1 
Pi /Pi = -2.28 ± 0.94 Myr"^ for Star 2 

For both Cepheids P/P is determined only with 2a accu- 
racy. Therefore, we do not consider the above values to be 
significant detections, but rather to provide observational 
limits on P. The results, nevertheless, are very interest- 
ing, because these limits seem to be in conflict with the 
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model predictions. The measured P/P values differ from 
theoretical ones (Tables 4 and 5) by 3.8f7 for Star 1 and 
by 4.6(7 for Star 2, respectively. If real, such a discrepancy 
would mean serious troubles for the evolutionary models. 
However, we should treat this result with a grain of salt. 
The quoted errors of P/P arc formal statistical errors re- 
sulting from the nonlinear least square procedure. Such 
errors very often underestimate the true uncertainities of 
the fited parameters. This is especially true in the case of 
multiperiodic pulsations (Costa, Kepler & Winget 1999; 
Costa & Kepler 2000). 

The P determinations reported in this section require 
confirmation with better data. This will soon be possi- 
ble when OGLE III photometry is released. Combined 
OGLE II+III datasct will provide the timcbase similar 
to that of MACHO, but much higher data quality should 
result in P errors about twice smaller. 

Further accumulation of data will eventually allow de- 
terminaton of P for all three modes. This is quite essential 
for understanding the true nature of the period changes. If 
they are caused by stellar evolution, then all modes should 
yield roughly the same values of P/P. No such equality 
is expected if period changes result from nonlinear mode 
interactions. 

7. Conclusions and discussion 

We show that pulsation periods measured in two LMC 
triple-mode Cepheids are consistent with values calculated 
for the first three overtones only if the envelope metal 
abundance is between Z = 0.004 and 0.007. This is some- 
what lower than Z = 0.008, typically adopted for the 
LMC. In one of the stars, fitting the periods requires the 
envelope hydrogen abundance X ^ 0.72, which is again 
somewhat below the commonly accepted LMC value of 
X = 0.74. The conclusions regarding envelope chemical 
composition rests mainly on the measured period ratios. 
For calculating the model periods, we rely on the linear 
pulsation theory. Wc believe, that nonlinear corrections 
to the periods arc small, but we cannot estimate them. 
Therefore, it is important that the derived metallicity lim- 
its are verified through direct spectroscopic determination. 

Effective temperatures implied by the models are con- 
sistent with observed {V — I) colours of the two stars. 
Unfortunately, the accuracy of the latter is currently too 
low to yield useful! constraints. Our models have lumi- 
nosities confined to narrow ranges. The inferred distance 
moduh are 18^"34 - 18™43 for one star and 18™49 - 18™53 
for the other. The difference between the two objects can 
be in large part attributed to photometric errors. The de- 
rived distances to the triple-mode Cepheids are within the 
range currently considered for the LMC, but favour rather 
longer LMC distances in this range. 

We find that parameters inferred from the pulsation 
periods of the two stars are consistent with evolutionary 
models undergoing the Hertzsprung gap crossing. The con- 
sistency does not require any overshooting from the con- 
vective core during the main sequence evolution. Seismic 



models do not exclude overshooting, though, but rather 
impose an upper limit on its extent. The limit is partic- 
ularly strong in case of LMC SC5 338399. The extent of 
overshooting can be further constrained once metal abun- 
dances or, alternatively, mean deredenned colours of the 
stars are presisely measured. 

Both the periods and the luminosities of the triple- 
mode Cepheids cannot be reconciled with evolutionary 
models in the phase of central helium burning. There 
is nothing strange in our proposal that these stars are 
undergoing Hertzsprung gap crossing. Stars with masses 
of ~ 3Mq, which we find for our objects, spend about 
10^ years crossing the instability strip on their way to the 
red giant branch. This is comparable to time spent inside 
the strip by more massive stars (M ~ 7Mq) during the 
central helium burning phase. The pulsation periods of 
such stars are between 10 and 20 days. We observe many 
LMC Cepheids in this period range. 

What may be a problem for the first crossing scenario 
is a possible conflict between predicted and measured rates 
of period change. The currently available data are insuf- 
ficient for significant P determination. We note, however, 
that in both Cepheids the most probable P values are 
more than 3.8a away from theoretical expectations. Since 
we are somewhat sceptical about P measurement errors, 
we do not yet consider this result to invalidate our models. 
But this is a warning sign. 

If the future determinations confirm the discrepancy, 
our models will be in serious troubles. In our opinion, the 
problem would be more severe than the one encountered 
by Pietrukowicz (2003) for long period Cepheids. He has 
found that theoretically predicted rates of period change 
are faster than the measured ones, but the size of the effect 
depends on whose evolutionary tracks are used. According 
to our unpublished calculations, there is no significant con- 
fiict if one uses the tracks of Girardi et al. (2000). 

The speed of evolution during the second and third 
crossings of the instability strip is affected by the element 
mixing in three distinct evolutionary phases. There are 
uncertainties regarding mixing efficiency in each of them. 
During the main sequence phase, we do not know the 
extent of the convectivc overshooting from the hydrogen 
burning core. During the red giant phase, we do not know 
the efficiency of the dredge up of the helium enriched ma- 
terial into the deep convective envelope. Finally, during 
the blue loop phase, there is a subtle problem with treat- 
ment of the boundary of the helium burning core. In the 
c;ase of the; first crossing of the instability strip, only the 
uncertainty of the main sequence overshooting affects the 
results, and this has been taken into account in our cal- 
culations. We believe that our prediction for the rates of 
period change is quite robust. Therefore, the possible dis- 
crepancy between calculated and observed P would pose a 
very serious challenge to our understanding of stellar evo- 
lution. Further observations of both triple-mode Cepheids 
are needed to settle this important issue. 
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